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ABSTRACT : We consider a real-time multi-server System with identical servers (such as machine controllers, 
unmanned aerial vehicles,overhearing devices, etc.) which can be adjusted/programmed for different types of 
activities (e.g. active or passive). This system provides a service for real-time jobs arriving via several channels 
(such as assembly lines, surveillance regions, communication ‘channels, etc.) and involves maintenance. We 
perform the worst case analysis of the system \vorking under maximum load with preemptive priorities assigned 
for servers of different activity type. We consider a system with separate queue to each channel. Two models 
with ample maintenance teams and shortage of maintenance teams are treated. We provide analytical 
approximations of steady state probabilities for these real-time systems and check their quality. 
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4 .-I. Introduction 1 ’— 2 

Real-time systems (RTS) are imbedded in modem technological structures, such as iproduction control 
systems, telecommunications systems, self-guided missiles, radars, £ ircraft, reconnaissance, etc. 

. «p_i 7?-. - *• 

According to Stankovic [1]: “real time systems are thosesystems that produce results in a timely 
manner'”, i.e. an action performed out of time constraints (too early or too late) may be useless, or harmful - 
even if such an action or computation is functionally correct. 


\ 


Input and output of RTS are mostly uncertain or incomplete. Many of RTS include therefore stochastic as 
well as dynamic components. —"V* 

' ri, l 1 - 

' r’ \ \ • i i i i I 

There exists a rich literature covering various real time models. 

We will focus on RTS with a zero deadline for Hie beginning of job processing. In these systems, jobs are 
processed immediately upon arrival, conditional on system availability. That part of the job which is not 
processed immediately is lost forever, and queueing of jobs (or. their parts) in these systems is impossible. 

The following works study this type of RTS. Kreimer and Mehrez ([2],‘ [3] ) proved that the non-mix policy 
maximizes the availability of a multiserver single-channel RTS involving preventive maintenance and working 
in general regime with any arrival pattern under consideration and .constant service and maintenance times. In 
[4] and [5] multi-server and multi-channel RTS (with ample and, 1 limited number of maintenance facilities 
respectively), working under maximum load regime were treated as finite source queues ([6]). In [7] various 
performance measures for RTS with arbitrary number of channels operating under a maximum load regime 
were studied. In [8£ [9^ <and [10] multi-server -and-multi-channel RTS working in general regime were 
presented. 


In [11] we have shown that even very large number of servers in RTS with ample maintenance facilities does 
not guarantee the maximum system availability, and optimal routing probabilities were computed analytically 
(for exponentially distributed service times) and via Cross Entropy (CE) [12],[13], [14] simulation approach 
(for generally distributed service times). These results were extended for RTS with limited maintenance 
facilities in [15]. 

RTS with priorities were studied in [16], [17], [18] (preemptive) and [19] (non-preemptive) correspondingly. 
In [16] a multi-server and multi-channel RTS with separate queues of servers for each channel and preemptive 
priorities was studied and a set of balance linear equations for steady-state probabilities was obtained. 
Unfortunately, these equations do not have analytical product-form solutions. 

In [19] several approximations methods, using the modifications of techniques proposed in [20], [21]were 
tested, and the best one was chosen. 
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The work presented here provides analytical approximations (based on the best method proposed in [18]) of 
steady state probabilities for RTS with separate queues of servers for each channel and preemptive priorities. 
We compare the approximation results with exact values of steady state probabilities and system performances. 

The paper is organized as follows: In Section 2, the description of the model with preemptive priorities is 
presented. In Section 3 balance equations (see [16]) for models with ample maintenance teams and shortage of 
maintenance teams are presented. Section 4 provides analytical approximations for these models. In Section 5 
some numerical results are presented. Finally, Section 6 is devoted to conclusions. 


II. Description of the model 

The most important characteristics of RTS with a zero deadline for the beginning of job processing are 
summarized in [8]. A real-world problem was studied in [16]. Here we provide the formal description of RTS 
from [16]. 

The system consists of r identical channels. For proper performance each channel needs exactly one fixed 
server at any moment (worst case), otherwise the information in this channel (at this specific moment) is lost. 

There are N servers (which are subject to breakdowns) in the system. A server, which is out of order, needs If 
time units of maintenance. After repair a fixed server may be of u-th type of quality(u=l,...,m) and is assigned 
to the v-th channel with probability p uv , (u=l, ...,m;v=l, These probabilities can be used as control 
parameters. Only after the repair is completed, the quality control procedure determines the quality type of 
fixed server. The fixed server of«-th type assigned to v-th channel is operative for a period of time S u before 

requiring Rj hours of repair. S u and If areindependent exponentially distributed random values with 

parameters p (u=l,...,m; v=I,...,r) and A respectively. It is assumed that there are A" identical maintenance 
teams in the system. Each team can repair exactly one server at a time. We will consider two models: 1) with 
ample maintenance facilities ( K > N ) so that all N servers can be repaired simultaneously if necessary, and 2) 
with shortage of maintenance facilities (K<N), in which case some of broken servers will wait for maintenance. 

The duration times R j of repair are i.r.v. exponentially distributed with parameter A , which does not 
depend on the quality type of the server (neither before nor after the repair). After repair, the fixed server will 
either be on stand-by or operating inside the channel. There is a separate queue of servers to each channel. 

We assume that servers of the first kind of quality type have the highest priority, servers of the second 
quality type are the next priority in line, and so on. Finally, servers of the m-th quality type have the lowest 
priority. Server operating in any channel is interrupted , if another fixed server of higher priority type arrives 
from maintenance. When the operating server must be repaired, the fixed server with highest priority takes its 
place. 

The system works under a maximum load (worst case) of nonstop data arrival to each channel). Thus, there 
exist a total of exactly r jobs in the whole system at any moment, and the nonstop operation of each channel is 
needed. 

If, during some period of time of length T .there is no fixed server to provide the proper operation of the 
channel, we will say that the part of the job of length T is lost forever. 


III. Balance equations for steady state probabilities 


In [15] the state of the system was defined as 


f ^ 


n lr ,...,n m r 


, where n u v u=l,...,m; v = 1, ...,r is a number 


m r 

of fixed servers of M-th activity type assigned for v-th channel (obviously, ^ ^ n < N and n > 0), and 

M=lv=l 
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P \ the corresponding steady state probability. There are 

"1,11 


f N + m-r 


states in total. We also denote 


yi,r’-"' n m,r ; 

?. u v = the rate of assignments of fixed servers of w-th activity type to the v-th channel. 

We will consider two models: 1) with ample maintenance facilities (K > N ); and 2) with shortage of 
maintenance facilities ( K<N). 


We will focus on the case r < N , otherwise all fixed servers will be busy and preemptive priority regime, 
therefore, will not work. 


3.1 Model with ample maintenance facilities 

We assume that there are ample identical maintenance facilities K > N available to repair alW servers 
simultaneously, if needed. Thus, each broken server enter maintenance facilities without delay. This model was 
studied, and the following set of linear equations for steady state probabilities was obtained in [16]. 


P f \ 


n l,v’—’ n u,v’".’ n m,v 
V U \ ,r’-"’ n m,r ) 

r m 


f r m ^ 

r m 

f 

( M Y 

o,i-^Z»i,v 

\ 


a+XZ 1 ™ 11 

n t v , max 

Mt,v 

/ 

_V .=! ,=i J 

v=lf=l 


^ i= 1 J, 


r m / x 

= ZZmin(n (i>v ,l)-P r 


v=lw=l 


+ min 


c „ „ ^ 

n \ ,l>—> n m ,1 


n.i v ,...,n uv l,...,n mv 

^ n l,r’”-’ n m,r J 
\ r m 


Af + l-I 2>,., 

v=li=l J 


+ 


-Hi} 

v v=li=l y v=\u=\ 


f n „ \ -min 

n l, l’ — ’ n m ,1 


”l,v. n u, v+h->”m,v 

n L,r’—’ n m,r 


n u v + 1, max 


u -1 


0,1-$, In, 


M« 


*'=i J) 


( 1 ) 


r N N-n h i N 

Z Z Z ... Z p, 

V=1«1,1=0 "2,1=0 n m = 0 


= 1 , 


" 1 , 1 . 
«l,v 
n l,r> • 


( 2 ) 


r m [0,f ^1] [0,W^1| 

where A = ZZ4,v^r = 1 “=1. . . NZ = 1 “J^t =\. , 

v=l«=l l,t = l 1, M = 1 


3.2 Model with shortage of maintenance facilities 

Here we assume that there are K (K < N ) maintenance facilities in the system. Thus, a shortage of 
maintenance facilities is possible when there are more than K broken servers in the system. In that case, the 
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broken server waits for maintenance. This model was studied in [16], and the following set of linear balance 
equations for steady state probabilities was obtained. 


p 

t n u> . t 


n \ v ,...,n uv ,...,n mv 


V n l,r’—’ n m,r y 



r r m ^ 

r m 

r 

f m ^ 

\ 

min 

K,N-^Tn itV ,K 

^+ZZ min 

n tv , max 

o,i-^z««> 

Mt, v 


K v=li=l ) 

v=n=t 

V 

^ i=l ), 

/ 


= XZmin(n UiV ,l )•/> 


v=1m=1 


/ >. min 


«i,v>•••>««, n m , v 

v n l 


^ r m "\ 

v v=t/=t y 


4 ; v + 



r r m N 

r m 



r 


+ min 


■ZZ p 


•min 

n uv +l,max 

oj-^Z^v 


V=1 i=l y 

v=lu=l 



V 

v i=i jj 


”l,v >•••>"»,v+h->”m,v 
n \,r ’-•’ n m,r 


Mu 


(3) 


r N N ~ n u N 

Z Z Z - Z p 

V=1»1,1=0 n 2 , 1 =0 n m =0 


«1,V 

v «l,r 


>”m,l 
-«m,v 
■> n m,r 


= 1(4) 


IV. Approximations 

We apply the Method 3 developed in [18] to get analytical approximations for solutions of equations (l)-(2). 
This method is a combination of the methods presented in [20] and [21]. These approximations instead of 
Global Balance Equations ((l)-(2) or (3)-(4)) use Local Balance Equations, and provide analytical product 
form solutions (see [6]) for steady state probabilities. We will use an RSS (root of sum of squares) criteria in 
order to evaluate the quality of these approximations: 


m -1 

N-Z'Tv 

r N N ~r>u ,=i 

RSS = SQRT{Y J Z Z Z [Pf 


”1,1.«m,l 

w l,v V5 w m,v 


-P\ ^] 2 }. 

n l,v ,...,n mv 
n\ r ,...,n, r 


y n l,r’-’ n m,r ^ \ n l,r’-’ n m,v 


(5) 


where P, 


"l.l.-.'Vl 


y n l,r’—’ n m,r 


and P 


r, l,l >—>«m,l 


\Jl\,r ’—> n m,r 

analytical approximation respectively. 


are exact numerical solution (of equations (l)-(2) or (3)-(4)) and 


Then analytical approximations of steady state probabilities for the model with ample maintenance teams and 
separate queues of servers to each channel are as follows. 
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f( r m \ \ 
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The corresponding analytical approximations of steady state probabilities for the model with shortage of 
maintenance teams and separate queues of servers to each channel are as follows. 
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where n, = max 


mm 


1 =1 i=m 


1 

o 

; S u „ = < 

hn u<v =0] 

y 

u,v 

0,n UfV >0 ] 


and 


m -1 

r N H—n 

EE E ■■■ If 

v=l«u = 0«2,i=0 n m ,v =0 


n l,v v 

n \ 


= 1, from which the value of PV n n can be easily found 


fo.0 

0.0 


V. Numerical results 

In this Section, we present some numerical results, showing the quality of our approximations in terms of RSS. 
Table 1 contains exact values of steady state probabilities for the model with ample maintenance facilities 
obtained from equations (l)-(2) presented in [16]. 

We consider the model with following parameters r = 2; IV = 3; 1 = 28; m = 2; u,, = 8, u?i = 10, Un = 6, an = 
^',Pn = 0.21,p 2 i = 0.29, P i 2 = 0.14,^22 = 0.36; K>N ’ 


Table 1. Steady state probabilities - exact numerical solution of equations (l)-(2) 


njj,n 2 i 

ni2,n 22 

P(n) 

iin,n 2 i 

n l2 , n 22 

P(n) 

ni 2 ,n 22 

P(n) 

nn,n 21 

n 12 ,n 22 

P(n) 

0,0 

0.008 

0,0 

0.003 

1,0 

0.003 

1,1 

0.045 

0,0 


3,0 


2,0 


0,1 


0,0 

0.027 

1,0 

0.009 

0,1 

0.044 

1,1 

0.021 

0,1 


0,0 


0,2 


1,0 


0,0 

0.005 

0,1 

0.028 

0,1 

0.038 

0,2 

0.046 

1,0 


0,0 


1,1 


0,1 


0,0 

0.058 

1,0 

0.024 

0,1 

0.007 

0,2 

0.024 

0,2 


0,1 


2,0 


1,0 


0,0 

0.029 

1,0 

0.005 

2,0 

0.009 

3,0 

0.007 

1,1 


1,0 


0,0 


0,0 


0,0 

0.004 

0,1 

0.051 

1,1 

0.039 

2,1 

0.034 

2,0 


0,1 


0,0 


0,0 


0,0 

0.048 

0,1 

0.017 

0,2 

0.060 

1,2 

0.084 

0,3 


1,0 


0,0 


0,0 


0,0 

0.087 

1,0 

0.030 

2,0 

0.012 

0,3 

0.048 

1,2 


0,2 


0,1 


0,0 


0,0 

0.027 

1,0 

0.016 

2,0 

0.005 



2,1 


1,1 


1,0 





Table 2 contains their analytical approximations counterparts for the same parameters: r = 2; /V = 3; /' = 28; m = 
2; nu = 8, fi 2 i = 10, fi 12 = 6, /t 2 2 = 12; p u = 0.21, p 2I = 0.29, p 12 = 0.14, p 22 = 0.36; K>N. 
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Table 2 .Steady state probabilities - analytical approximations. 


nu,n 2 i 

11/2,1122 

P'(n) 

11/1,1121 

n !2, n 22 

nil) 

11/1,112/ 

11/2,1122 

P'(n) 

11/1,112/ 

11/2,1122 

P'(n) 

0,0 

0.016 

0,0 

0.005 

1,0 

0.004 

U 

0.021 

0,0 


3,0 


2,0 


0,1 


0,0 

0.040 

1,0 

0.012 

0,1 

0.053 

1,1 

0.017 

0,1 


0,0 


0,2 


1,0 


0,0 

0.007 

0,1 

0.038 

0,1 

0.043 

0,2 

0.051 

1,0 


0,0 


1,1 


0,1 


0,0 

0.067 

1,0 

0.020 

0,1 

0.013 

0,2 

0.041 

0,2 


0,1 


2,0 


1,0 


0,0 

0.020 

1,0 

0.006 

2,0 

0.009 

3,0 

0.007 

1,1 


1,0 


0,0 


0,0 


0,0 

0.004 

0,1 

0.064 

1,1 

0.026 

2,1 

0.019 

2,0 


0,1 


0,0 


0,0 


0,0 

0.056 

0,1 

0.019 

0,2 

0.061 

1,2 

0.046 

0,3 


1,0 


0,0 


0,0 


0,0 

0.111 

1,0 

0.017 

2,0 

0.012 

0,3 

0.049 

1,2 


0,2 


0,1 


0,0 


0,0 

0.013 

1,0 

0.013 

2,0 

0.004 



2,1 


1,1 


1,0 





Corresponding RSSE=0.068 (between results of Table 1 and Table 2) 


Table 3 contains exact numerical values as well as analytical approximations of several Performance 
characteristics (averages) for the same model: r = 2; N = 3; 1 = 28; m = 2; / i u = 8, /i 2 ; = 10, /in = 6, jin = 12; pu 
= 0.21, p 2I = 0.29, p 12 = 0. 14, p 22 = 0.36; K>N. 


Table 3 . Exact values and analytical approximations of Performance characteristics (averages) 


Performance characteristic 

Exact 

Value 

Analytical 

Approximations 

Nu. of fixed servers L 

2.542 

2.461 

Nu. of fixed serversat channel 1 /V, | } 

1.308 

1.142 

Nu. of fixed serversat channel 2 /Vp ; 

1.233 

1.319 

Nu. of fixed serversof type 1 N ] 

0.748 

0.646 

Nu. of fixed serversof type 2 N 

1.793 

1.815 

Nu. of fixed serversof type 1 at channel 1 j + Q 11 

0.412 

0.282 

Nu. of fixed serversof type 2 at channel 1 L ? t + Q , j 

0.896 

0.860 

Nu. of fixed serversof type 1 at channel 2 Lj 0 + <2i -> 

0.336 

0.364 
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Nu. of fixed serversof type 2 at channel 2L 1? + Q 1 2 

0.897 

0.955 

Availability at channel 1 AVj 

0.725 

0.661 

Availability at channel 2 A\> 2 

0.650 

0.696 

Availability of the system Av 

0.688 

0.688 


VI. Conclusions 

We have found a good analytical approximations for RTS with a separate queues for each channel, and 
preemptive priorities (for servers of different quality type) working under maximal load. Two models with 
ample maintenance teams and shortage of maintenance teams were treated and numerical results showing the 
quality of approximations are provided. 
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